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Abstract. Let .4 be a class of objects, equipped with an integer size such that for all n the 
number a„ of objects of size n is finite. We are interested in the case where the generating 
function X^^^ a„f" is rational, or more generally algebraic. This property has a practical 
interest, since one can usually say a lot on the numbers a^, but also a combinatorial one: 
the rational or algebraic nature of the generating function suggests that the objects have 
a (possibly hidden) structure, similar to the linear structure of words in the rational case, 
and to the branching structure of trees in the algebraic case. We describe and illustrate 
this combinatorial intuition, and discuss its validity. While it seems to be satisfactory in 
the rational case, it is probably incomplete in the algebraic one. We conclude with open 
questions. 
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1. Introduction 

The general topic of this paper is the enumeration of discrete objects (words, trees, 
graphs...) and more specifically the rational or algebraic nature of the associated 
generating functions. Let ^ be a class of discrete objects equipped with a size: 

size : A ^ 

A ^ \A\. 

Assume that for all n, the number a„ of objects of size n is finite. The generating 
function of the objects of A, counted by their size, is the following formal power 
series in the indeterminate t: 

A(t):=^a„i"=^il^l. (1) 

n>0 AgA 

To take a very simple example, if A is the set of words on the alphabet {a, b} 
and the size of a word is its number of letters, then the generating function is 

E„>o2"i" = l/(l-20. 

Generating functions provide both a tool for solving counting problems, and a 
concise way to encode their solution. Ideally, one would probably dream of finding 
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a closed formula for the numbers a„. But the world of mathematical objects would 
be extremely poor if this was always possible. In practise, one is usually happy with 
an expression of the generating function A{t), or even with a recurrence relation 
defining the sequence a„, or a functional equation defining A{t). 

Enumerative problems arise spontaneously in various fields of mathematics, 
computer science, and physics. Among the most generous suppliers of such prob- 
lems, let us cite discrete probability theory, the analysis of the complexity of algo- 
rithms |56[l44j. and the discrete models of statistical physics, like the famous Ising 
model [5]. More generally, counting the objects that occur in one's work seems 
to answer a natural curiosity. It helps to understand the objects, for instance to 
appreciate how restrictive are the conditions that define them. It also forces us 
to get some understanding of the structure of the objects: an enumerative result 
never comes for free, but only after one has elucidated, at least partly, what the 
objects really are. 

We focus in this survey on objects having a rational, or, more generally, al- 
gebraic generating function. Rational and algebraic formal power series are well- 
behaved objects with many interesting properties. This is one of the reasons why 
several classical textbooks on enumeration devote one or several chapters to these 
series [43l[7H[75]. These chapters give typical examples of objects with a rational 
[resp. algebraic] generating function (GF). After a while, the collection of these 
examples builds up a general picture: one starts thinking that yes, all these ob- 
jects have something in common in their structure. At the same time arises the 
following question: do all objects with a rational [algebraic] GF look like that? In 
other words, what does it mean, what does it suggest about the objects when they 
are counted by a rational [algebraic] GF ? 

This question is at the heart of this survey. For each of the two classes of series 
under consideration, we first present a general family of enumerative problems 
whose solution falls invariably in this class. These problems are simple to describe: 
the first one deals with walks in a directed graph, the other with plane trees. 
Interestingly, these families of objects admit alternative descriptions in language 
theoretic terms: they correspond to regular languages, and to unambiguous context- 
free languages, respectively. The words of these languages have a clear recursive 
structure, which explains directly the rationality [algebraicity] of their GF. 

The series counting words of a regular [unambiguous context-free] language are 
called N-rational [N-algebraic]. It is worth noting that a rational [algebraic] series 
with non-negative coefficients is not necessarily N-rational [N-algebraic]. Since we 
want to appreciate whether our two generic classes of objects are good represen- 
tatives of objects with a rational [algebraic] GF, the first question to address is 
the following: do we always fall in the class of N-rational [N-algebraic] series when 
we count objects with a rational [algebraic] GF? More informally, do these objects 
exhibit a structure similar to the structure of regular [context-free] languages? Is 
such a structure usually clearly visible? That is to say, is it easy to feel, to predict 
rationality [algebraicity]? 

We shall see that the answer to all these questions tends to be yes in the rational 
case (with a few warnings...) but is probably no in the algebraic case. In particular. 
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the rich world of planar maps (planar graphs embedded in the sphere) abounds in 
candidates for non-N-algebraicity. The algebraicity of the associated GFs has been 
known for more than 40 years (at least for some families of maps) , but it is only 
in the past 10 years that a general combinatorial explanation of this algebraicity 
has emerged. Moreover, the underlying constructions are more general that those 
allowed in context-free descriptions, as they involve taking complements. 

Each of the main two sections ends with a Hst of questions. In particular, we 
present at the end of Section [3] several counting problems that are simple to state 
and have an algebraic GF, but for reasons that remain mysterious. 

The paper is sometimes written in an informal style. We hope that this will 
not stop the reader. We have tried to give precise references where he/she will find 
more details and more material on the topics we discuss. In particular, this survey 
borrows a lot to two books that we warmly recommend: Stanley's Enumerative 
Combinatorics \74\ 175). and Flajolet & Sedgewick's Analytic Combinatorics |43) . 

Notation and definitions. Given a (commutative) ring R, we denote by R[t] 
the ring of polynomials in t having coefficients in R. A Laurent series in t is a 
series of the form A{t) = J2n>no '^"t" , with uq G Z and an & R for all n. If 
uq > 0, we say that A{t) is a formal power series. The coefficient of t" is denoted 
a„ := [t"]yl(i). The set of Laurent series forms a ring, and even a field if i? is a 
field. The quasi-inverse of A{t) is the series A*{t) := 1/(1 — A{t)). If A{t) is a 
formal power series with constant term 0, then A*{t) is a formal power series too. 

In most occasions, the series we consider are GFs of the form ^ and thus 
have rational coefficients. However, we sometimes consider refined enumeration 
problems, in which every object A is weighted, usually by a monomial w{A) in some 
additional indeterminates xi, . . . ,Xm- The weighted GF is then X^Ae^ '"^(^)^'"^' ' 
so that the coefficient ring is Q[xi, . . . , Xm] rather than Q. 

We denote Jfc| = {1, 2, . . . , A:}. We use the standard notation N, Z, Q, and 
P:= {1,2,3,...}. 



2. Rational generating functions 

2.1. Definitions and properties. The Laurent series A{t) with coeffi- 
cients in the field R is said to be rational if it can be written in the form 

A(t) ^ 
^ ^ Q{t) 

where P{t) and Q{t) belong to R[t]. 

There is probably no need to spend a lot of time explaining why such series are 
simple and well-behaved. We refer to [74^ Ch. 4] and |431 Ch. IV] for a survey of 
their properties. Let us review briefiy some of them, in the case where R = Q. The 
set of (Laurent) rational series is closed under sum, product, derivation, reciprocals 
- but not under integration as shown by A{t) = 1/(1 — i). The coefficients a„ of 
a rational series A{t) satisfy a linear recurrence relation with constant coefficients: 
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for n large enough, 

a„ = cia„_i + C2a„_2 H h Cfea„_fc. 

The partial fraction expansion of A{t) provides a closed form expression of these 
coefficients of the form: 

k 

a„ = ^P,;(nK (2) 

where the /ii are the reciprocals of the roots of the denominator Q{t), and the Pi are 
polynomials. In particular, if A{t) has non-negative integer coefficients, its radius 
of convergence p is one its the poles (Pringsheim) and the "typical" asymptotic 
behaviour of a„ is 

a„ - np-'^n'^ (3) 

where G N and k is an algebraic number. The above statement has to be taken 
with a grain of salt: all poles of minimal modulus may actually contribute to the 
dominant term in the asymptotic expansion of a„ , as indicated by ^ . 

Let us add that Fade approximants allow us to guess whether a generating 
function whose first coefficients are known is Hkely to be rational. For instance, 
given the 10 first coefficients of the series 

A{t) = t+2 i^+e 1®+6466 1^+20727 1^°+0 (t") 

it is easy to conjecture that actually 



A{t) 



1 - 5 1 + 7 <2 - 4 <3 ■ 



Fade approximants are implemented in most computer algebra packages. For in- 
stance, the relevant Maple command is convert/ratpoly. 

2.2. Walks on a digraph. We now introduce our typical "rational" objects. 
Let G — (y, E) be a directed graph with (finite) vertex set V = \p\ and (directed) 
edge set E <Z V X V . A walk of length n on G is a sequence of vertices w = 
(iiQ, f 1, . . . , Vn) such that for all i, the pair {vi, Wi+i) is an edge. Such a walk goes 
from vq to Vn. We denote \w\ = n. Now assign to each directed edge e a weight 
(an indeterminate) Xe- Define the weight Xw of the walk w as the product of the 
weights of the edges it visits: more precisely, 

1=0 

See Fig. 1(a) for an example. Let X denote the (weighted) adjacency matrix of G: 
for i and j in the entry Xij is Xe if («, j) = e is an edge of G and otherwise. 
Let Wij{t) be the weighted generating function of walks going from i to j: 

W^,J{t) = X^t^^l 
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Figure 1. (a) A weighted digraph. The default value of the weight is 1. (b) A deterministic 
automaton on the alphabet {a, b, c, d, b, c}. The initial state is 1 and and the final states 
are 2 and 3. 



It is well-known, and easy to prove, that Wij is a rational function in t with 
coefEcicients in Q[xe,e G E] (see O Thm. 4.7.1]). 

Theorem 2.1. The series Wij{t) is the [i, i)-entry in the matrix (1 — tX)'"^ . 

This theorem reduces the enumeration of walks on a digraph to the calculation 
of the inverse of a matrix with polynomial coefficients. It seems to be little known in 
the combinatorics community that this inverse matrix can be computed by studying 
the elementary cycles of the digraph G. This practical tool relies on Viennot's 
theory of heaps of pieces [81]. Since it is little known, and often convenient, let us 
advertise it here. It will be illustrated further down. 

An elementary cycle of G is a closed walk w ~ {vq,vi, . . . , f„_i, vq) such that 
Vf), . . . , Vn-i are distinct. It is defined up to a cyclic permutation of the Vi. That 
is, {vi,V2, . . . , Vn-i,vo, vi) is the same cycle as w. A collection 7 = {71, . . . , 7, } of 
(elementary) cycles is non-intersecting if the 7^ are pairwise disjoint. The weight 
Xj of 7 is the product of the weights of the 7^. We denote I7I = l7i|- 

Proposition 2.2 ([H]). The generating function of walks going from i to j reads 
where 

D= J2 i-^Vx^t^^^ and N^^j = ^ {~lYx^x^t\^\+\^\. 

7 = {7i,...,7r} w;'i = {ji,...,jr} 

The polynomial D is the alternating generating function of non-intersecting collec- 
tions of cycles. In the expression of N , j a non-intersecting collection of cycles 
and w a self-avoiding walk going from i to j , disjoint from the cycles of 



To illustrate this result, let us determine the generating function of walks going 
from 1 to 2 and from 1 to 3 on the digraph of Fig. [DJa). This graph contains 4 
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cycles of length 1, 2 cycles of length 2, 2 cycles of length 3 and 1 cycle of length 4. 
By forming all non-intersecting collections of cycles, one finds: 

D{t) = l-{3 + x)t + {3 + 3x-2)f + {-l-3x + 3 + x-2)t^ + {x-l-x + l + x-y)t'^ 
= 1 - (3 + x)i + (1 + 3x)t'^ - 2xt^ + {x - y)t^. 

There is only one self-avoiding walk (SAW) going from 1 to 2, and one SAW going 
from 1 to 3 (via the vertex 2). The collections of cycles that do not intersect these 
walks are formed of loops, which gives 

Ni^2 = til-tf{l-xt) and Ni^a = t^{l - tf. 

Hence the generating function of walks that start from 1 and end at 2 or 3 is: 

_ A^i,2 + ^i.3 ^ tji-tra+t-xt) 

"^1,2-^^^1,3 l^{3 + x)t+{l + 3x)t^-2xt^ + {x-y)t'^' ^' 

2.3. Regular languages and automata. There is a very close connec- 
tion between the collection of walks on a digraph and the words of regular lan- 
guages. Let A be an alphabet, that is, a finite set of symbols (called letters). A 
word on ^ is a sequence u = uiU2 ■ ■ ■ Un of letters. The number of occurrences of 
the letter a in the word u is denoted \u\a. The product of two words uiU2 ■ ■ - Un and 
V1V2 • • • is the concatenation U1U2 • • • UnViV2 • ■ -Vm- The empty word is denoted 
e. A language on ^ is a set of words. We define two operations on languages: 

- the product CIC of two languages £ and IC is the set of words uv, with u G C 
and V g IC; this product is easily seen to be associative, 

- the star C* of the language C is the union of all languages , for fc > 0. By 
convention, £° is reduced to the empty word e. 

A finite state automaton on ^ is a digraph (V, E) with possibly multiple edges, 
together with: 

- a labelling of the edges by letters of A, that is to say, a function L : E ^ A, 

- an initial vertex i, 

- a set Vf C F of final vertices. 

The vertices are usually called the states of the automaton. The automaton is 
deterministic if for every state v and every letter a, there is at most one edge 
labelled a starting from v. 

To every walk on the underlying multigraph, one associates a word on the 
alphabet A by reading the letters met along the walk. The language C recognized 
by the automaton is the set of words associated with walks going from the initial 
state i to one of the states of Vf. For j € V, let Cj denote the set of words 
associated with walks going from i to j. These sets admit a recursive description. 
For the automaton of Fig. [Ifb), one has £ = £2 U £3 with 

A = {e}, 

£2 = £ic U £20 U £30 U £40, £4 = £20 U £30 U £40 U £56, 

£3 = £2CU£35U£3C, £5 = £2cU£35u£3cU£56. 
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Remarkably, there also exists a non-recursive combinatorial description of the 
languages that are recognized by an automaton [52l Thms. 3.3 and 3.10]. 

Theorem 2.3. Let C be a language on the alphabet A. There exists a finite state 
automaton that recognizes C if and only if C can be expressed in terms of finite 
languages on A, using a finite number of unions, products and stars of languages. 

If these conditions hold, C is said to be regular. Moreover, there exists a de- 
terministic automaton that recognizes C. 

Regular languages and walks on digraphs. Take a deterministic automaton, 
and associate with it a weighted digraph as follows: the vertices are those of the 
automaton, and for all vertices j and fc, if m edges go from j to k in the automaton, 
they are replaced by a single edge labelled m in the digraph. For instance, the 
automaton of Fig. [T]Jb) gives the digraph to its left, with x = y = 2. Clearly, the 
length GF of words of C is the GF of (weighted) walks of this digraph going from 
the initial vertex i to one of the final vertices of Vf. For instance, according to 
the length GF of the language recognized by the automaton of Fig. [ijb) is 



Take a regular language £ recognized by a deterministic automaton A. There 
exists another deterministic automaton that recognizes C and does not contain 
multiple edges. The key is to create a state {j, a) for every edge labelled a ending 
at j in the automaton A. The digraph associated with this new automaton has all 
its edges labelled 1, so that there exists a length preserving bijection between the 
words of £ and the walks on the digraph going from a specified initial vertex vq to 
one of the vertices of a given subset Vf of vertices. 

Conversely, starting from a digraph with all edges labelled 1, together with a 
specified vertex vq and a set Vf of final vertices, it is easy to construct a regular 
language that is in bijection with the walks of the graph going from vq to Vf 
(consider the automaton obtained by labeUing all edges with distinct letters). This 
shows that counting words of regular languages is completely equivalent to counting 
walks in digraphs. In particular, the set of rational series obtained in both types 
of problems coincide, and have even been given a name: 

Definition 2.4. A series A{t) ~ J2n>o "with coefficients in N is said to be 
N-rational if there exists a regular language having generating function A{t) — oq. 

The description of regular languages given by Theorem 12 .31 implies that the set 
of N-rational series contains the smallest set of series containing N[t] and closed 
under sum, product and quasi-inverse. The converse is true [HI Thm. II. 5.1]. 
There exists a simple way to decide whether a given rational series with coefficients 
in N is N-rational [TTl Thms. II.10.2 and II.10.5]. 

Theorem 2.5. A series A{t) = Sri>o ^^'^ coefficients in N is N-rational if 
and only if there exists a positive integer p such that for all r e {0, . . . ,p}, the 
series 



A{t) 



t{l-tf 



(5) 



l-5t + 7t2-4t3' 




n>0 
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has a unique singularity of minimal modulus (called dominant^. 

There exist rational series with non-negative integer coefficients that are not 
N-rational. For instance, let a be such that cos a = 3/5 and sin a = 4/5, and 
define a„ ~ 25"cos(na)^. It is not hard to see that a„ is a non-negative integer. 
The associated series A{t) reads 

^ l-2t + 225t^ 

^' (1 - 25<)(625<2 + 14t+ !)■ 

It has 3 distinct dominant poles. As a is not a rational multiple of tt, the same 
holds for all series Ao^p{t), for all values of p. Thus A{t) is not N-rational. 

2.4. The combinatorial intuition of rational generating func- 
tions. We have described two families of combinatorial objects that naturally 
yield rational generating functions: walks in a digraph and words of regular lan- 
guages. We have, moreover, shown that the enumeration of these objects are 
equivalent problems. It seems that these families convey the "right" intuition about 
objects with a rational GF. By this, we mean informally that: 

(i) "every" family of objects with a rational GF has actually an N-rational GF, 

(m) for almost all families of combinatorial objects with a rational GF, it is easy 
to foresee that there will be a bijection between these objects and words of 
a regular language. 

Point (ii) means that most of these families J- have a clear automatic structure, 
similar to the automatic structure of regular languages: roughly speaking, the 
objects of T can be constructed recursively using unions of sets and concatenation 
of cells (replacing letters). A more formal definition would simply paraphrase the 
definition of automata. 

Point (i) means simply that I have never met a counting problem that would 
yield a rational, but not N-rational GF. This includes problems coming from al- 
gebra, like growth functions of groups. On the contrary. Point (n) only concerns 
purely combinatorial problems (but I do not want to be asked about the border 
between combinatorics and algebra). It admits very few counter-examples. Some 
will be discussed in Section [231 For the moment, let us illustrate the two above 
statements by describing the automatic structure of certain classes of objects (some 
being rather general), borrowed from [74l Ch. 4]. 

2.4.1. Column-convex polyominoes. A polyomino is a finite union of cells 
of the square lattice, whose interior is connected. Polyominoes are considered 
up to a translation. A polyomino is column-convex (cc) if its intersection with 
every vertical line is connected. Let a„ be the number of cc-polyominoes having 
n cells, and let A(t) be the associated generating function. We claim that these 
polyominoes have an automatic structure. 

Consider a cc-polyomino P having n cells. Let us number these cells from 1 to 
n as illustrated in Fig. [51 The columns are visited from left to right. In the first 
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Figure 2. A column-convex polyomino, with the numbering and encoding of the cells. 



column, cells are numbered from bottom to top. In each of the other columns, the 
lowest cell that has a left neighbour gets the smallest number; then the cells lying 
below it are numbered from top to bottom, and finally the cells lying above it are 
numbered from bottom to top. Note that for all i, the cells labelled 1, 2, . . . , i form 
a cc-polyomino. This will be essential in our description of the automatic structure 
of these objects. Associate with P the word u = mi • • • w„ on the alphabet {a, b, c} 
defined by 

- Ui = c (like Column) if the zth cell is the first to be visited in its column, 

- Ui = b (like Below) if the ith cell Hes below the first visited cell of its column, 
-Ui = a (like Above) if the ith cell lies above the first visited cell of its column. 

Then, add a bar on the letter Ui if the ith cell of P has a South neighbour, an 
East neighbour, but no South-East neighbour. (In other words, the barred letters 
indicate where to start a new column, when the bottommost cell of this new 
column Hes above the bottommost cell of the previous column.) This gives a word 
V on the alphabet {a,b,c,a,b,c}. It is not hard to see that the map that sends 
P on the word w is a size-preserving bijection between cc-polyominoes and words 
recognized by the automaton of Fig. [Ijb) . Hence by ((5]) , the generating function 
of column-convex polyominoes is [76] : 

_ t{l-tf 
^ ' ~ l-5< + 7t2-4t3- 

2.4.2. P-partitions. A partition of the integer n into at most k parts is a non- 
decreasing fc-tuple A = (Ai, . . . , Afe) of nonnegative integers that sum to n. This 
classical number-theoretic notion is generaHzed by the notion of P-partitions. Let 
P be a natural partial order on \k\ (by natural we mean that if i < j in P, then 
z < j in N). A P -partition of ti is a fc-tuple A = (Ai, . . . , Afc) of nonnegative integers 
that sum to n and satisfy Ai < Xj ii i < j in P. Thus when P is the natural total 
order on a P-partition is simply a partitioifl. 
We are interested in the following series: 

Pp(t) = ^tl^ 

A 

where the sum runs over all P-partitions and |A| = Ai + • • • + A^ is the weight of A. 



P-partition is usually defined as an order-reversing map from [fc| to N |741 Section 4.5]. 
Both notions are of course completely equivalent. 
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The case of ordinary partitions is easy to analyze: every partition can be written 
in a unique way as a linear combination 

ciA(i) + --- + CfeA« (6) 

where A^*^ = (0, 0, . . . , 0, 1, 1, . . . , 1) has exactly i parts equal to 1 and a G N. The 
weight of A^*-' is i, and one obtains: 

^^^^^ ^ (i-t){i-t^)---{i-t^y 

The automatic structure of (ordinary) partitions is transparent: since they are 
constructed by adding a number of copies of A^^-*, then a number of copies of A'^^^ 
and so on, there is a size preserving bijection between these partitions and walks 
starting from 1 and ending anywhere in the following digraph: 




[4] 

Note that this graph corresponds to fc = 4, and that an edge labelled [£] must be 
understood as a sequence of £ edges. These labels do not correspond to multipHci- 
ties. Observe that the only cycles in this digraph are loops. This, combined with 
Proposition 12.21 explains the factored form of the denominator of ([7]) . 

Consider now the partial order on |4] defined by 1 < 3, 2 < 3 and 2 < 4. The 
partitions of weight at most 2 are 



(0,0,0,0), (0,0, 1,0), (0,0,0,1), (1,0, 1,0), (0,0, 1,1), (0,0, 2,0), (0,0, 0,2), 

so that Fp{t) = 1 + 2t + At^ + 0{t^) . If one is brave enough to list P-partitions of 
weight at most 20, the Fade approximant of the truncated series thus obtained is 
remarkably simple: 

(l-i)(l-t2)(l-i3)(l_t4)+^(^ 

and allows one to make a (correct) conjecture. 

It turns out that the generating function of P-partitions is always a rational 
series of denominator (1 — t){l — . . . {1 — t^). Moreover, P-partitions obey our 
general intuition about objects with a rational GF. The following proposition, 
illustrated below by an example, describes their automatic structure: the set of 
P-partitions can be partitioned into a finite number of subsets; in each of these 
subsets, partitions have a structure similar to ([6]). Recall that a linear extension 
of P is a bijection a on \k\ such that < a{j) if i < j in P. 
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Proposition 2.6 ([74], Section 4.5). Let P he a natural order on \k\. 

For every P -partition A, there exists a unique linear extension o of P such that 
for all i, Xa{i) < Acr(i+i), the inequality being strict if a{i) > a{i + 1). We say that 
A is compatible with a . 

Given a linear extension a, the P-partitions that are compatible with a can be 
written in a unique way as a linear combination with coefficients in N; 

A('^'")+ciA('^^i) + --- + CfeA('^^'=' (8) 

where A'-'^'^-' is the smallest P -partition compatible with a: 

and for 1 < i < k, 

has exactly i parts equal to 1. Thus the GF of these P-partitions is 

^e(cr) 

^^'"^^^ " (l-t)(l-i2)...(l_tfc) 

where e(a) is a variant of the Major index of a: 

e(a)= (^-*)- 

i:CT(i)>(T(i+l) 

Example. Let us return to the order 1 < 3, 2 < 3 and 2 < 4. The 5 linear 
extensions are 1234, 2134, 1243, 2143 and 2413. Take cr = 2143. The P-partitions A 
that are compatible with a are those that satisfy A2 < Ai < A4 < A3. The smallest 
of them is thus A^'"-") = (1,0,2,1). Then A^""'^' = (0,0,1,0), A^'^'^) ^ (0,0,1,1), 
^ 0, 1, 1) and A('^'4) ^ (1, 1, 1, 1). 

2.4.3. Integer points in a convex polyhedral cone (|74|. Sec. 4.6). Let H 

be a finite collection of Hnear half-spaces of M™ of the form aai + • • • + Cma„i > 0, 
with a G Z. We are interested in the set £ of non-negative integer points a ~ 
(ai, . . . , am) lying in the intersection of those half-spaces. For instance, we could 
have the following set £, illustrated in Fig. [31(a): 

£ = {(ai, 0:2) e : 2ai > 02 and 2q:2 > aj. (9) 

Numerous enumerative problems (including P-partitions) can be formulated in 
terms of linear inequalities as above. The generating function of £ is 

i?(t) = ^tH, 

where \a\ = ai -\ ha™. In the above example, E{t) = 1 + 1'^ + 2t^ + 1"^ + 2t^ + 

St^ + 2t'^ + 0{t^). 
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The set £ is a monoid (it is closed under summation). In general, it is not a free 
monoid. Geometrically, the set C of non-negative real points in the intersection 
of the half-spaces of H forms a pointed convex polyhedral cone (the term pointed 
means that it does not contain a line), and £ is the set of integer points in C. 

The simplicial case. In the simplest case, the cone C is simplicial. This implies 
that the monoid £ is simplicial, meaning that there exists linearly independent 
vectors a^^\ . . . , a'-'"'^ such that 



a 



qia 



(1) 



gfca^^' with G Q, > o} 



This is the case in Example ([9]), with a^^^ = (1, 2) and a^^^ = (2, 1). The interior 
of £ (the set of points of £ that are not on the boundary of C) is then 



£ 



a e N™ : a = q^a^^' + ■ ■ ■ + qt 



a^^'^ with qi e Q, q, > o| 



(10) 



Then there exists a finite subset V of £ [resp. V of £] such that every element of 
£ [resp. £] can be written uniquely in the form 



/3 + Cia(l) + hcfco; 



(k) 



(11) 



with (3 £ V [resp. (3 £ V] and c; G N |741 Lemma 4.6.7]. In our running exam- 
ple 111) Jtaken with a^^) = (1, 2) and a^^) = (2, 1), one has V = {(0, 0), (1, 1), (2, 2)} 
while V = {(1, 1), (2,2), (3J3)}. Compare ^ with the structure found for P- 
partitions ^ . Thus £ and £ have an automatic structure and their GFs read 



E{t)^ 



nti(i-ti"''i) 



resp. E{t) 



nti(i-ii""'i) 



In Example (|9]), one thus obtains 

Eit) ' + + ' 



t + t' 



(l-<3)2 (l-t)(l-i3) 



and E{t) = t^E{t). 
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The general case. The set £ can always be partitioned into a finite number of 
sets T of the form |[10]), where T \& a. simpUcial monoid [TU Ch. 4, Eq. (24)]. Thus 
£, as a finite union of sets with an automatic structure, has an automatic structure 
as well. The associated generating function E{t) is N-rational, with a denominator 
which is a product of cyclotomic polynomials. 
Consider, for example, the set 

E = {(«!, q;2, as) e : ^3 < ai + Q!2}. 

The cone C of non-negative real points a satisfying <ai+ is not simphcial, 
as it has 4 faces of dimension 2, lying respectively in the hyperplanes ai — for 
1 = 1,2,3 and as = ai + 0.2 (Fig. [31(b)). But it is the union of two simplicial cones 
C\ and C2, obtained by intersecting C with the half-spaces a\ > 03 and ai < a^, 
respectively. Let £1 [resp. £2] denote the set of integer points of Ci [resp. C2]. 
The fastest way to obtain the generating function E{t) is to write 

E{t) = Eiit) + E2it) - Ei2{t) (12) 

where Ei2{t) counts integer points in the intersection of Ci and C2 (that is, in 
the plane ai = a^). Since £1, £2 and £1 n £2 are simplicial cones (of dimension 
3, 3 and 2 respectively), the method presented above for simplicial cones applies. 
Indeed, £1 [resp. £2; £12] is the set of hnear combinations (with coefficients in N) 
of (1,0, 1), (0,1,0) and (1,0,0) [resp. (1,0,1), (0,1,0) and (0,1,1); (1,0,1) and 
(0,1,0)]. This implies: 



E{t) 



1 l + t + t^ 



(l-i)2(l-f2) (l_^)(l_i2)2 (l-f)(l-<2) (l-t)(l-f2)2 



However, the "minus" sign in (fT2| prevents us from seeing directly the automatic 
nature of £ (the difference of N-rational series is not always N-rational). This 
structure only becomes clear when we write £ as the disjoint union of the interiors 
of all simplicial monoids induced by the triangulation of C into Ci and C2. These 
monoids are the integer points of the faces (of all possible dimensions) of Ci and 
C2. As there are 12 such faces (more precisely, 1 [resp. 4, 5, 2] faces of dimension 
[resp. 1, 2, 3]), this gives £ as the disjoint union of 12 sets having an automatic 
structure of the form (fTOl) . 



2.5. Rational generating functions: more difficult questions. 



2.5.1. Predicting rationality. We wrote in Section [2^ that it is usually easy to 
foresee, to predict when a class of combinatorial objects has a rational GF. There 
are a few exceptions. Here is one of the most remarl-cable ones. 

Example 2.7 (Directed animals). A directed animal with a compact source of 
size A: is a finite set of points A on the square lattice 1? such that: 

- the points (— i, i) for < z < k belong to A; they are called the source points, 

- all the other points in A can be reached from one of the source points by a 
path made of North and East steps, having all its vertices in A. 
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Figure 4. Compact-source directed animals on the square and triangular lattices. 

See Fig. |4] for an illustration. A similar notion exists for the triangular lattice. It 
turns out that these animals have extremely simple generating functions [50^ llOj. 

Theorem 2.8. The number of compact-source directed animals of cardinality n is 
3""^ on the square lattice, and 4"~^ on the triangular lattice. 

The corresponding GFs are respectively t/(l — 3t) and — At), and are 
as rational as a series can be. There is at the moment no simple combinatorial 
intuition as to why these animals have rational GFs. A bijection between square 
lattice animals and words on a 3- letter alphabet was described in [50], but it does 
not shed a clear light on the structure of these objects. Still, there is now a 
convincing explanation of the algebraicity of these series (see Section 13.4. 2p . 

Example 2.9 (The area under Dyck paths). Another family of (slightly less 
natural) examples is provided by the enumeration of points lying below certain 
lattice paths. For instance, let us call Dyck path of length 2n any path P on 
formed of steps (1, 1) and (1, —1), that starts from (0, 0) and ends at (2n, 0) without 
ever hitting a point with a negative ordinate. The area below P is the number of 
non-negative integer points (i, j), with i < 2n, lying weakly below P (Fig. [5|). It 
turns out that the sum of the areas of Dyck paths of length 2n is simply 

P:\P\=2n 

Again, the rationality of the associated generating function does not seem easy to 
predict, but there are good combinatorial reasons explaining why it is algebraic. 
See [33t I65j for a direct explanation of this result, references, and a few variations 
on this phenomenon, first spotted by Kreweras [58] . 




Figure 5. The 5 Dyck paths of length 6 and the 4"^ — 64 points lying below. 



Finally, let us mention that our optimistic statement about how easy it is to 
predict the rationality of a generating function becomes less and less true as we 
move from purely combinatorial problems to more algebraic ones. For instance, it 
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is not especially easy to foresee that a group has an automatic structure [39]. Let 
us give also an example coming from number theory. Let P{x) = P{xi, . . . , Xr) be 
a polynomial with integer coefficients, and take p a prime. For n > 0, let a„ be 
the number of a; € (Z/p"Z)'' such that P{x) = modp". Then the generating 
function X)n '^n^" is rational. A related result holds with p-adic solutions [371153] . 

2.5.2. Computing a rational generating function. Let us start with an ele- 
mentary, but important observation. Many enumerative problems, including some 
very hard, can be approximated by problems having a rational GF. To take one ex- 
ample, consider the notoriously difficult problem of counting self-avoiding polygons 
(elementary cycles) on the square lattice. It is easy to convince oneself that the 
generating function of SAP lying in a horizontal strip of height k is rational for all 
k. This does not mean that it will be easy (or even possible, in the current state of 
affairs) to compute the corresponding generating function when k = 100. Needless 
to say, there is at the moment no hope to express this GF for a generic value of k. 
The generating function of SAP having 2k horizontal steps can also be seen to be 
rational. Moreover, these SAP can be described in terms of Hnear inequalities (as 
in Section I2.4.3P , which implies that the denominator of the corresponding series 
Gk is a product of cyclotomic polynomials. But again, no one knows what this 
series is for a generic value of k, or even for k = 100. Still, some progress have been 
made recently, since it has been proved that the series Gk have more and more 
poles as k increases, which means that their denominators involve infinitely many 
cyclotomic polynomials [68] . This may be considered as a proof of the difficulty of 
this enumerative problem [51j . 

In general, computing the (rational) generating function of a family of objects 
depending on a parameter k may be non-obvious, if not difficult, even if the objects 
are clearly regular, and even if the final result turns out to be nice. A classical 
example is provided by the growth functions of Coxeter groups [61]. Here is a more 
combinatorial example. A partition A = (Ai, . . . , Afe) is said to be a k-Lecture Hall 
partition (fc-LHP) if 

Ai Ao Xk 
0< — < — <•••<— . 
- 1 - 2 - - k 

Since these partitions are defined by linear inequalities, it follows from Section [2.4.3l 

that their weight generating function is rational, with a denominator formed of 

cyclotomic polynomials. Still, there is no clear reason to expect that |15j : 



Several proofs have been given for this result and variations on it. See for in- 
stance [16t [35] and references in the latter paper. Some of these proofs are based 
on a bijection between lecture hall partitions and partitions into parts taken in 
{1, 3, . . . , 2fc — 1}, but these bijections are never really simple [82j |4Q]. 

2.5.3. N-rationality. As we wrote in Section [531 we do not know of a counting 
problem that would yield a rational, but not N-rational series. It would certainly 
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be interesting to find one (even if it ruins some parts of this paper) . 

Let us return to Soittola's criterion for N-rationality (Theorem 12. Sp . It is not 
always easy to prove that a rational series has non-negative coefficients. For in- 
stance, it was conjectured in [4^ that for any odd fc, the number of partitions of 
n into parts taken in {fc, A: -I- 1, . . . , 2fc — 1} is a non-decreasing function of n, for 
n > 1. In terms of generating functions, this means that the series 

, i_q 

(l-(zfe)(l-(7fe+i)---(l-<?2fe-i) 

has non-negative coefficients. This was only proved recently [67]. When k is even, 
a similar result holds for the series 

i__q 

^ (1 - - '7'''+^) • • • (1 - q^'')i^ - q^^^^) ■ 

Once the non-negativity of the coefficients has been established, it is not hard to 
prove that these series are N-rational. This raises the question of finding a family 
of combinatorial objects that they count. 



3. Algebraic generating functions 

3.1. Definitions and properties. The Laurent series A{t) with coeffi- 
cients in the field R is said to be algebraic (over R{t)) if it satisfies a non-trivial 
algebraic equation: 

P{t,A{t)) = Q 

where P is a bivariate polynomial with coefficients in R. 

We assume below that i? = Q. Again, the set of algebraic Laurent series 
possesses numerous interesting properties [75t Ch. 6], [43^ Ch. VII]. It is closed 
under sum, product, derivation, reciprocals, but not under integration. These clo- 
sure properties become effective using either the theory of elimination or Grobner 
bases, which are implemented in most computer algebra packages. The coefficients 
a„ of an algebraic series A{t) satisfy a Hnear recurrence relation with polynomial 
coefficients: for n large enough, 

Po(?^)ar^ -l-pi(n)a„_i +p2{n)an~2 H ^ Pk{n)an-k = 0. 

Thus the first n coefficients can be computed using a linear number of operations. 

There is no systematic way to express the coefficients of an algebraic series in 
closed form. Still, one can sometimes apply the Lagrange inversion formula: 

Proposition 3.1. Let <& and ^' be two formal power series and let U = U{t) be 
the unique formal power series with no constant term satisfying 

U = t^{U). 

Then for n > 0, the coefficient oft" in ^'(C/) is: 

[t"]^(C/) = i[t"-i](^'(i)$(t)"). 



Rational and algebraic series in combinatorial enumeration 



17 



Given an algebraic equation P(t, A{t)) = 0, one can decide whether there exists 
a series U{t) and two rational series $ and ^' satisfying 

U^t^iU) and A^^{U). (13) 

Indeed, such series exist if and only if the genus of the curve P{t, a) is zero [U 
Ch. 15]. Moreover, both the genus and a parametrization of the curve in the 
form flS]) can be determined algorithmically. 

Example 3.2 (Finding a rational parametrization). The following algebraic 
equation was recently obtained [22], after a highly non-combinatorial derivation, 
for the GF of certain planar graphs carrying a hard-particle configuration: 

= 23328 t^A^ + 27 1^ (91 - 2088 1) A^ + 1^ (86 - 3951 1 + 46710 + 3456 1^) A^ 

+ (l -mt+ 1598 - 11743 - 14544 1^) A-1+66 i-1495 1^+11485 1^+128 1'^. 

(14) 

The package algcurves of Maple, and more precisely the commands genus and 
parametrization, reveal that a rational parametrization is obtained by setting 

(3 + 7) (9 [/2 + 33 C/ + 37) 

ij - — - 3 i . 

{■iu + if 

Of course, this is just the net result of Maple, which is not necessarily very 
meaningful for combinatorics. Still, starting from this parametrization, one obtains 
after a few attempts an alternative parametrizing series V with positive coefficients: 



(l-2y)(l-3y + 3y2)- ^ 

The main interest of such a parametrization for this problem does not lie in the 
possibility of applying the Lagrange inversion formula. Rather, it suggests that a 
more combinatorial approach exists, based on the enumeration of certain trees, in 
the vein of [191 [27]. It also gives a hint of what these trees may look like. 

Another convenient tool borrowed from the theory of algebraic curves is the 
possibility to explore all branches of the curve P(t, A{t)) = in the neighbourhood 
of a given point t^. This is based on Newton's polygon method. All branches have 
a Puiseux expansion, that is, an expansion of the form: 

Ait)^ ^ a„(i-io)"/' 

n>no 

with no G Z, fi € P. The coefficients a„ belong to C (in general, to an algebraic clo- 
sure of the ground field) . These expansions can be computed automatically using 
standard software. For instance, the Maple command puiseux of the algcurves 
package tells us that lfT4|) has a unique solution that is a formal power series, the 
other three solutions starting with a term i"^. 
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Such Puiseux expansions are crucial for studying the asymptotic behaviour of 
the coefRcients of an algebraic series A{t). As in the rational case, one has first 
to locate the singularities of A{t), considered as a function of a complex variable 
t. These singularities are found among the roots of the discriminant and of the 
leading coefficient of P{t, a) (seen as a polynomial in a). The singular expansion of 
A{t) near its singularities of smallest modulus can then be converted, using certain 
transfer theorems, into an asymptotic expansion of the coefficients [42 | l43 t VII. 4]. 

Example 3.3 (Asymptotics of the coefRcients of an algebraic series). 

Consider the series V{t) defined by (fTSl) . Its singularities lie among the roots of 
the discriminant 

A{t) = -3 + 114t - 4635t^ + 55296t^ 

Only one root is real. Denote it to ^ 0.065. The modulus of the other two roots is 
smaller than to, so they could, in theory, be candidates for singularities. However, 
V{t) has non-negative coefficients, and this implies, by Pringsheim's theorem, that 
one of the roots of minimal modulus is real and positive. Hence V{t) has a unique 
singularity, lying at to- A Puiseux expansion at this point gives 



V{t) = Co - civ/1 - t/to + 0{t - to), 
for some explicit (positive) algebraic numbers cq and ci, which translates into 

r]V{t) = ^t~"n-'/'il + o{l)). 

The determination of asymptotic expansions for the coefficients of algebraic 
series is probably not far from being completely automated, at least in the case of 
series with non-negative coefficients [3ll . The "typical" behaviour is 

P"""', (16) 



r(fi + i) 



where k is an algebraic number and e Q \ { — 1, —2, —3, . . .}. Compare with the 
result ^ obtained for rational series. Again, the above statement is not exact, as 
the contribution of all dominant singularities must be taken into account. See O 
Thm. VII. 6] for a complete statement. 

Let us add that, again, one can guess if a series A{t) given by its first coefficients 
satisfies an algebraic equation P{t,A{t)) = of a given bi-degree {d,e). The 
guessing procedure requires to know at least {d+l){e + l) coefficients, and amounts 
to solving a system of linear equations. It is implemented in the package Gf un of 
Maple [T^- For instance, given the 10 first coefficients of the series V(t) satisfying 
V{0) = and lfT5|) . one automatically conjectures ifTS]). 

3.2. Plane trees. Our typical "algebraic" objects will be (plane) trees. Let 
us begin with their usual intuitive recursive definition. A tree is a graph formed 
of a distinguished vertex (called the root) to which are attached a certain number 
(possibly zero) of trees, ordered from left to right. The number of these trees is 
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Figure 6. (a) A plane tree, (b) A rooted planar map. (c) The corresponding 4-valent 
map (thick lines). 



the degree of the root. The roots of these trees are the children of the root. A 
more rigorous definition describes a tree as a finite set of words on the alphabet 
P satisfying certain conditions [63]. We hope that our less formal definition and 
Fig. [6l[a) suffice to understand what we mean. The vertices of a tree are often 
called nodes. Nodes of degree are called leaves, the others are called inner nodes. 

The enumeration of classes of trees yields very often algebraic equations. Let 
us consider for instance the complete binary trees, that is, the trees in which all 
vertices have degree or 2 (Fig. [I2l) . Let an be the number of such trees having n 
leaves. Then, by looking at the two (sub)trees of the root, one gets, for n > 1: 

n-l 

fe=l 

The initial condition is ai = 1. In terms of GFs, this gives A{t) = t + A{t)^, which 
is easily solved: 



Ait) 



2 ^ n+l 

Tl>0 



(17) 



More generally, many algebraic series obtained in enumeration are given as the 
first component of the solution of a system of the form 



A, = P,{t,Ai,...,Ak), 



(18) 



for some polynomials Pi{t,xi, . . . ,Xk) having coefficients in Z. This system is said 
to be proper if Pi has no constant term [Pi{Q, ■ • ■ , 0) =0) and does not contain 
any linear term Xi. It is positive if the coefficients of the Pi are non-negative. For 
instance, 

Ai=t^ + A1A2 and A2 = 2tAl 

is a proper positive system. The system is quadratic if every Pi{t,xi, . . . , Xk) is a 
linear combination of the monomials t and xgXm, for 1 < £ < m < k. 

Theorem 3.4 ([75], Thm. 6.6.10 and [71], Thm. IV.2.2). A proper algebraic 
system has a unique solution (Ai, . . . ,Ak) in the set of formal power series in t 
with no constant term. This solution is called the canonical solution of the system. 
The series Ai is also the first component of the solution of 
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- a proper quadratic system, 

- a proper system of the form Bi = tQi{t, Bi, . . . , Bi), for 1 < i < i. 

These two systems can be chosen to be positive if the original system is positive. 

Proof. Let us prove the last property, which we have not found in the above refer- 
ences. Assume Ai satisfies ifTSj) and that this system is quadratic. The ith equation 
reads Ai = rriit + ni^£r(j-)A^(j). Rewrite each monomial AiAj as tUij and add the 

equations J7y = t [mirrij + minjUcr(j)r{j) + mjn,Uc,(i)r{i) + '^i'^3Ua[i)T{i)Ua{j)T{j)) ■ 
The new system has the required properties. □ 

Definition 3.5. A series A(t) is N-algebraic if it has coefficients in N and if A{t) — 
A{Q) is the first component of the solution of a proper positive system. 

Proper positive systems like l|18p can always be given a combinatorial inter- 
pretation in terms of trees. Every vertex of these trees carries a label {i, c) where 
i G fk} and c g P. We say that i is the type of the vertex and that c is its colour. The 
type of a tree is the type of its root. Write Aq = t, so that Ai = ^^(^o, ^i, ■ • ■ , ^fc)- 
Let ^0 be the set reduced to the tree with one node, labelled (0, 1). For z G 
let Ai be the set of trees such that 

- the root has type i, 

- the types of the subtrees of the root, visited from left to right, are 0, . . . , 0, 1, . . . , 
1, . . . ,k, . . . ,k, in this order, 

- if exactly Cj children of the root have type j, the colour of the root is any inte- 
ger in the interval [1, m], where m is the coefficient of Xq° ■ ■ ■ x'j.'' in Pi{xo, . . . ,Xk). 

Then it is not hard to see that Ai (t) is the generating function of trees of type 
i, counted by the number of leaves. This explains why trees will be, in the rest of 
this paper, our typical "algebraic" objects. 

3.3. Context-free languages. As in the case of rational (and, more pre- 
cisely, N-rational) series, there exists a family of languages that is closely related 
to algebraic series. A context-free grammar G consists of 

- a set 5 = {5*1, . . . , Sk} of symbols, with one distinguished symbol, say. Si, 

- a finite alphabet A of letters, disjoint from S, 

- a set of rewriting rules of the form Si ^ w where w is a non-empty word on 
the alphabet S U A. 

The grammar is proper if there is no rule Si — *■ Sj . The language £{G) generated 
by G is the set of words on the alphabet A that can be obtained from Si by applying 
iteratively the rewriting rules. A language is context-free is there exists a context- 
free grammar that generates it. In this case there exists also a proper context-free 
grammar that generates it. 

Example 3.6 (Dyek words). Consider the grammar G having only one symbol, 
S, alphabet {a, 6}, and rules S ^ ab + abS + aSb + aSbS (which is short for 
S ^ ab, 5" — *■ abS, S aSb, S aSbS). It is easy to see that C{G) is the 
set of non-empty words u on {a, 6} such that \u\a = \u\b and for every prefix v of 
u, \v\a > \v\b. These words, called Dyck words, provide a simple encoding of the 
Dyck paths met in Example 12.91 
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A derivation tree associated with G is a plane tree in which all inner nodes 
are labelled by symbols, and all leaves by letters, in such a way that if a node is 
labelled Si and its children wi,. . . ,Wk (from left to right), then the rewriting rule 
Si wi - ■ ■ Wk is in the grammar. If the root is labelled Si, then the word obtained 
by reading the labels of the leaves in prefix order {i.e., from left to right) belongs 
to the language generated by G. Conversely, for every word w in C{G), there exists 
at least one derivation tree with root labelled Si that gives w. The grammar is 
said to be unambiguous if every word of C{G) admits a unique derivation tree. 

Assume G is proper. For 1 < i < k, let Ai(t) be the generating function of 
derivation trees rooted at Si, counted by the number of leaves. With each rule r, 
associate the monomial M(r) = Xq" • • • x^*" where eo [resp. e^, with z > 0] is the 
number of letters of A [resp. occurrences of Si] in the right-hand side of r. Then 
the series Ai, . . . ,Ak form the canonical solution of the proper positive system (fTS)) . 
with 



where the sum runs over all rules r with left-hand side Si. 

Conversely, starting from a positive system Bi = tQi(t, Bi, . . . , Bk) and its 
canonical solution, it is always possible to construct an unambiguous grammar 
with symbols Si, . . . ,Sk such that Bi is the generating function of derivation trees 
rooted at Si (the idea is to introduce a new letter Ui for each occurrence of t). 
In view of Theorem 13.41 and Definition 13.51 this gives the following alternative 
characterization of N-algebraic series: 

Proposition 3.7. A series A{t) is N-algebraic if and only if only A{0) e N 
and there exists an unambiguous context-free language having generating function 



3.4. The combinatorial intuition of algebraic generating func- 
tions. We have described two families of combinatorial objects that naturally 
yield algebraic GFs: plane trees and words of unambiguous context-free languages. 
We have, moreover, shown a close relationship between these two types of objects. 
These two families convey the standard intuition of what a family with an algebraic 
generating function looks Hke: the algebraicity suggests that it may (or should...) 
be possible to give a recursive description of the objects based on disjoint union of 
sets and concatenation of objects. Underlying such a description is a context-free 
grammar. This intuition is the basis of the so-called Schiitzenberger methodology, 
according to which the "right" combinatorial way of proving algebraicity is to de- 
scribe a bijection between the objects one counts and the words of an unambiguous 
context-free language. This approach has led in the 80's and 90's to numerous sat- 
isfactory explanations of the algebraicity of certain series, and we describe some of 
them in this subsection. Let us, however, warn the reader that the similarities with 
the rational case will stop here. Indeed, it seems that the "context-free" intuition 
is far from explaining all algebraicity phenomena in enumerative combinatorics. In 
particular. 




A{t)~A{0). 
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(i) it is very likely that many families of objects have an algebraic, but not 
N-algebraic generating function, 

{ii) there are many families of combinatorial objects with an algebraic GF that 
do not exhibit a clear "context-free" structure, based on union and concate- 
nation. For several of these families, there is just no explanation of this type, 
be it clear or not. 

This will be discussed in the next subsections. For the moment, let us illustrate 
the "context-free" intuition. 

3.4.1. Walks on a line. Let 5 be a finite subset of Z. Let W be the set of walks 
on the fine Z that start from and take their steps in S. The length of a walk is its 
number of steps. Let Wk be the set of walks ending at position k. For k > 0, let 
A^fc be the subset of Wk consisting of walks that never visit a negative position, 
and let ^A be the union of the sets Mk- In probabilistic terms, the walks in M 
would be called meanders and the walks of ^Ao excursions. Of course, a walk is 
simply a sequence of steps, hence a word on the alphabet S. Thus the sets of walks 
we have defined can be considered as languages on this alphabet. 

Theorem 3.8. The language W is simply S* and is thus regular. The languages 
A4, Wk and Mk are unambiguous context-free for all k. 

Proof. We only describe the (very simple) case S = {+1,-1}, to illustrate the 
ideas that are involved in the construction of the grammar. We encode the steps 
by the letter a, the steps —1 by b, and introduce some auxiliary languages: 

• Mq , the subset of Wo formed of walks that never visit a positive position, 

• [resp. Wq], the subset of Wo formed of walks that start with a [resp. b]. 
The language Mo will be generated from the symbol Mq, and similarly for the 
other languages. By looking at the first time a walk of Mo [resp. Mq] reaches 
position after its first step, one obtains 

Mo ^ a{l + Mo)b{l + Mo) and Mf^ ^ b{l + Mf^)a{l + M^)- 

By considering the last visit to of a walk of A^fe, one obtains, for fc > 0: 

Mk ^ (l + Mo)a(U=i +A/fc_i). 

This is easily adapted to general meanders: 

A/^A/o + (l + il/o)a(l + Af). 

Considering the first step of a walk of Wo gives 

Wo + Wq with 1^0+ A/o(l + W^) and M^il + W^). 

Finally, for fc > 0, looking at the first visit at 1 [resp. —1] of a walk of Wk [resp. 
W-k] yields 

Wk^ il + M^)a{lk=i+Wk-i) [resp. W-k ^ (1 + Mo)b {lk=i + W^^k-i))] ■ 
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For a general set of steps S, various grammars have been described for the 
languages ^Ak oi meanders |38 [ [60 l [59], For Wk, we refer to |591 Section 4] where 
the (representative) case S = {—2,-1, 0, 1, 2} is treated. □ 



Theorem 13.81 is often described in terms of walks in starting from (0, 0) and 
taking their steps in G S}. The conditions on the positions of the walks 

that lead to the definition oiMk and Wk are restated in terms of conditions on the 
ordinates of the vertices visited by the walk. A harmless generahzation is obtained 
by taking steps in a finite subset 5 of P x Z. A walk is still encoded by a word 
on the alphabet S. The languages Wk remain unambiguous context-free. If each 
step is, moreover, weighted by a rational number Wij, then the generating 

function of walks of W, counted by the coordinates of their endpoint, is 

W{t,s)= ^ 



The generating function Wk{t) that counts (weighted) walks ending at ordinate 
k is the coefficient of in W{t,s). Since Wk is unambiguous context-free, the 
series Wk{t) is algebraic. This gives a combinatorial explanation of the fofiowing 
result [ZSl Thm. 6.3.3]. 

Theorem 3.9 (Diagonals of rational series). LetA{x,y) — X^m n>o '^™.n'^™2/" 
be a series in two variables x and y, with coefficients in Q, that is rational. Then 
the diagonal of A, that is, the series AA{t) = J2n>Q'^n.nt"' , is algebraic. 

Proof. By finearity, it suffices to consider the case 

At \ """y" 

A{x,y) 



1 J2o<rn,n<d ^m,nX"^y'" 

with co,o = 0. Set x = ts and y = t/s. The diagonal of A satisfies 

1 



AAit"^) = [s°]A(ts,t/s) t''+''[s''-''] 



1 _ „ fm+n „n 

^ l^O<m,n<d^™''"'' * 



which is algebraic as it counts weighted paths in Wb-a, for a certain set of steps. 
Hence AA{t) is algebraic too. □ 



The converse of Theorem [XQ] holds: every series B{t) that is algebraic over Q{t) 
is the diagonal of a bivariate rational series A{x,t) [70] . 

Note. If one is simply interested in obtaining a set of algebraic equations defining 
the GFs of the sets Mk and Wk, a more straightforward approach is to use a 
partial fraction decomposition (for Wk) and the kernel method (for Mk)- See [75l 
6.3], and [m Example 3]. 
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3.4.2. Directed animals. Let us move to an example where a neat context-free 
exists, but is uneasy to discover. We return to the directed animals defined in 
Section 12.5.11 As discussed there, there is no simple explanation as to why the 
number of compact-source animals is so simple (Theorem 12. Sp . Still, there is a 
convincing explanation for the algebraicity of the corresponding series: directed 
animals have, indeed, a context-free structure. This structure was discovered a 
few years after the proof of Theorem 12.81 with the development by Viennot of 
the theory of heaps [H], a geometric version of partially commutative monoids 
[30] . Intuitively, a heap is obtained by dropping vertically some sohd pieces, the 
one after the other. Thus, a piece lies either on the "floor" (then it is said to be 
minimal), or covers, at least partially, another piece. 

Directed animals are, in essence, heaps. To see this, replace every point of 
the animal by a dimer (Fig. [7]). Note that if the animal has a unique source, the 
associated heap has a unique minimal piece. Such heaps are named pyramids. 




Figure 7. (a) A directed animal and the associated pyramid, (b) A half-pyramid. 



What makes heaps interesting here is that there exists a monoid structure on 
the set of heaps: The product of two heaps is obtained by putting one heap above 
the other and dropping its pieces. This product is the key in our context-free 
description of directed animals. 

Let us begin with the description of pyramids (one-source animals) . A pyramid 
is either a half-pyramid (Fig. [lljb)), or the product of a half-pyramid and a pyramid 
(Fig. m top). Let P{t) denote the GF of pyramids counted by the number of 
dimers, and H{t) denote the GF of half-pyramids. Then P{t) = H{t){l + Pit)). 
Now, a half-pyramid may be reduced to a single dimer. If it has several dimers, it 
is the product of a single dimer and of one or two half-pyramids (Fig. [8l bottom), 
which implies H{t) = t + tH{t) + tH'^(t). 

A trivial computation finally provides the GF of directed (single-source) animals: 

(.MeHit)^'-!^^^^^^^^^^ 




2t I 

The enumeration of compact-source directed animals is equivalent to the enumer- 
ation of heaps having a compact basis (the minimal dimers are adjacent). The 
generating function of heaps having a compact basis formed with k dimers is 
P{t)H{t)^~^ (Fig. E]), which implies that the generating function of compact- 
source animals is 

P{t) t 
l~H{t) l-3t' 
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Figure 8. Decomposition of pyramids (top) and half-pyramids (bottom). 




Figure 9. Decomposition of heaps having a compact basis. 



3.5. The world of planar maps. We have seen in Section [321 that plane 
trees are the paradigm for objects with an algebraic generating function. A more 
general family of plane objects seems to be just as deeply associated with algebraic 
series, but for reasons that are far more mysterious: planar maps. 

A (planar) map is a proper embedding of a planar graph in the sphere (Fig.[6tb)). 
In order to avoid symmetries, all the maps we consider are rooted: this means that 
one edge is distinguished and oriented. Maps are only considered up to a continu- 
ous deformation of the sphere. A map induces a 2-cell decomposition of the sphere: 
the cells of dimension [resp. 1, 2] are called vertices [resp. edges, faces]. Hence 
plane trees are maps with a single face. 

The interest for the enumeration of planar maps dates back to the early 60 's, in 
connection with the 4-colour theorem. The first results are due to Tutte [TT j lTS l lTQ]. 
Ten to fifteen years later, maps started to be investigated independently in theoret- 
ical physics, as a model for 2-dimensional quantum gravity [28119]. However, neither 
the recursive approach used by Tutte and his disciples, nor the physics approach 
based on matrix integrals were able to explain in a combinatorially satisfactory 
way the following observations: 

- the generating functions of many classes of planar maps are algebraic, 

- the associated numbers are often irritatingly simple. 
Let us illustrate this with three examples: 

1. General maps. The number of planar maps having n edges is |80j : 




(19) 
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The associated generating function G — G{t) = X]n>o ffn^" satisfies: 

- 1 + 16t + (1 - m)G + 27t^G^ = 0. (20) 

2. Loopless triangulations. The number of loopless triangulations (maps in 
which all faces have degree 3) having 2n + 2 faces is [62] : 

2" /3n\ 
" (n + l)(2n+l) U A 

The associated generating function T = T{t) ~ J2n ^n^" satisfies 

1 - 27i + (-1 + 36<)T - 8tT'^ - 16fT^ = 0. 

3. Three-connected triangulations. The number of 3-connected triangulations 
having 2n + 2 faces is [77] : 

_ 2 Mn+ 1\ 

~ (n + l)(3n + 2) \ n )' 

The associated generating function M = M{t) = t„i" satisfies 

-1 + 16f + (1 - 2Qt)M + (3i + 8t^)A/2 + Si^M^ + t^Af'' = 0. 

These maps are in bijection with rooted maximal planar simple graphs (graphs 
with no loop nor multiple edge that lose planarity as soon as one adds an edge). 

At last, in the past ten years, a general combinatorial picture has emerged, 
suggesting that maps are, in essence, unrooted plane trees. In what follows, we 
illustrate on the example of general maps the main three approaches that now 
exist, and give references for further developments of these methods. 

3.5.1. The recursive approach. We leave to the reader to experience person- 
ally that maps do not have an obvious context-free structure. Still, maps do have 
a simple recursive structure, based on the deletion of the root-edge. However, in 
order to exploit this structure, one is forced to keep track of the degree of the root- 
face (the face lying to the right of the root edge) . The decomposition illustrated 
in Fig. [To] leads in a few lines to the following equation: 

GM ^l + tu'G{u,tr+tu ^G{u,t)-G{l,t) ^ 

u — \ 

where G{u, t) counts planar maps by the number of edges {t) and the degree of the 
root-face (u). 

It can be checked that the above equation defines G(u, t) uniquely as a formal 
power series in t (with polynomial coefficients in u). However, it is not clear on the 
equation why G{1, t) (and hence G(u, t)) are algebraic. In his original paper, Tutte 
first guessed the value of Gi(t) := G(l, t), and then proved the existence of a series 
G{u,t) that fits with Gi{t) when u = 1, and satisfies the above equation. Still, 
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+ 




Figure 10. Tutte's decomposition of rooted planar maps. 



a bit later, Brown came with a method for solving ((2T1) : the so-called quadratic 
method Sec. 2.9]. Write ^ in the form {2aG{u,t) + hf = 5, where a, 6 

and 5 are polynomials in t,u and Gi[t). That is, 

{2tu'^{u - l)G{u, t)+tu^ -u + lf ^ U^u^{u-l)Gi + {l~uf-Atu^+Qtu^+uH'^-2tu^. 

It is not hard to see, even without knowing the value of G(u,t), that there exists 
a (unique) formal power series in t, say U = U{t), that cancels the left-hand side 
of this equation. That is. 



This implies that the series C/ is a double root of the polynomial S that lies on the 
right-hand side. The discriminant of this polynomial (in u) thus vanishes: this 
gives the algebraic equation (|20|) satisfied by G(l,t). 

The enumeration of many other families of planar maps can also be attacked 
by a recursive description based on the deletion of an edge (or vertex, or face...). 
See for instance ^62j for 2-connected triangulations, or [6] for maps with prescribed 
face degrees. (For maps with high connectivity, like 3-connected triangulations, an 
additional composition formula is often required \77\ [3].) The resulting equations 
are usually of the form 



where F{u) = F{t, u), the main generating function, is a series in t with polynomial 
coefficients in u, and Fi,. . . ,Fk are series in t only, independent of u. Brown's 
quadratic method applies as long as the degree in F{u) is 2 (for the linear case, 
see the kernel method in [T7t[2|). Recently, it was understood how these equations 
could be solved in full generality [22]. Moreover, the solution of any (well-founded) 
equation of the above type was shown to be algebraic. This provides two types of 
enumerative results: 

- the proof that many map generating functions are algebraic: it now suffices 
to exhibit an equation of the form (|22| , or to explain why such an equation exists, 

- the solution of previously unsolved map problems (like the enumeration of 
hard-particle configurations on maps, which led to (fT4|) . or that of triangulations 
with high vertex degrees [8]). 

3.5.2. Matrix integrals. In the late 70's, it was understood by a group of physi- 
cists that certain matrix integral techniques coming from quantum field theory 
could be used to attack enumerative problems on maps pHl [9]. This approach 



U =l + tU^ + 2tU^{U - l)GiU, t). 



PiF{u),Fu...,Fk,t,u) = 0, 



(22) 
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proved to be extremely efficient (even if it is usually not fully rigorous) . The first 
step is fairly automatized, and consists in converting the description of maps into a 
certain integral. For instance, the relevant integral for the enumeration of 4-valent 
maps (maps in which all vertices have degree 4) is 

Z(t,N) = ^^-^J dHe y f / ), 

where the integration space is that of hermitian matrices H of size N, equipped 
with the Lebesgue measure dH = Y['^^kkY[k<e^^kedyke with hke = xm + iyki- 
As there is a classical bijection between 4-valent maps with n vertices and planar 
maps with n edges (Fig.[6Uc)), we are still dealing with our reference problem: the 
enumeration of general planax maps. The connection between the above integral 
and maps is 

G[t)^tE'(t) with E{t)= lim ■^logZit.N). 

N—*OQ N 

Other map problems lead to integrals involving several hermitian matrices [S^ . We 
refer to |83] for a neat explanation of the encoding of map problems by integrals, 
and to [45t [4T] (and references therein) for the evaluation of integrals. 

3.5.3. Planar maps and trees. We finally come to a combinatorial explanation 
of the formula/equation for Qn and G{t). Take a plane binary tree with n (inner) 
nodes, planted at a leaf, and add to every inner node a new distinguished child, 
called a hud. At each node, we have three choices for the position of the bud 
(Fig. [TlTa)). The new tree, called budding tree, has now n buds and n + 2 leaves. 
Now start from the root and walk around the tree in counterclockwise order, pay- 
ing attention to the sequence of buds and leaves you meet. Each time a bud is 
immediately followed by a leaf in this sequence, match them by forming a new edge 
(Fig.EUb)) and then go on walking around the plane figure thus obtained. At the 
end, exactly two leaves remain unmatched. Match them together and orient this 
final edge in one of the two possible ways. Also, mark the face to the left of the 
matching edge that ends at the root-leaf. 




Figure 11. (a). A budding tree, (b) An intermediate step in the matching procedure, (c) 
The resulting 4-valent map, with its marked face. 
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Theorem 3.10 ([73]). The above correspondence is a bijection between pairs (T, e) 
where T is a budding tree having n inner nodes and e E {0, 1}, and A-valent maps 
with n vertices and a marked face. 

The value of e tells how to orient the final matching edge. Schaeffer first used 
this bijection to explain combinatorially the formula ifig]). Indeed, the number of 
budding trees with n inner nodes is clearly 3"(^^")/(ri, + 1) (see l|17p ). while the 
number of 4-valent maps with n vertices and a marked face is [n + 2)gn. Eq. lfT9|) 
follows. 

Later, it was realized that this construction could also be used to explain the 
algebraicity of the series G{t) [23]. Say that a budding tree is balanced if the root- 
leaf is not matched by a bud. Take such a tree, match all buds, and orient the final 
edge from the root-leaf to the other unmatched leaf. This gives a bijection between 
balanced budding trees and 4-valent maps. We thus have to count balanced trees, 
or, equivalently, the unbalanced ones. By re-rooting them at the bud that matches 
the root-leaf, one sees that they are in bijection with a node attached to three 
budding trees. This gives 

G{t) = B{t) - tB{tf , where B{t) = it{l + B{t)f 

counts budding trees by (inner) nodes. The above construction involves taking 
a difference of N-algebraic series, which needs not be N-algebraic. We actually 
conjecture that the series G{t) is not N-algebraic (see Section [3. 6. 4p . 

There is little doubt that the above construction (once described in greater 
detail...) explains in a very satisfactory way both the simplicity of the formula 
giving gn and the algebraicity of G{t). Moreover, this is not an ad hoc, isolated 
magic trick: over the past ten years, it was realized that this construction is one in a 
family of constructions of the same type, which apply to numerous famihes of maps 
(Eulerian maps [73j, maps with prescribed vertex degrees [23], constellations |18j . 
bipartite maps with prescribed degrees [19], maps with higher connectivity |66| 
I47|). Definitely, these constructions reveal a lot about the combinatorial nature of 
planar maps. 

To conclude this section, let us mention that a different combinatorial construc- 
tion for general planar maps, discovered in the early 80's [34], has recently been 
simplified [32] and adapted to other famihes of maps [36l [Ml [25l [26] . It is a bit 
less easy to handle than the one based on trees with buds, but it allows one to 
keep track of the distances between some vertices of the map. This has led to re- 
markable connections with a random probability distribution called the Integrated 
SuperBrownian Excursion [32]. A third type of construction has emerged even 
more recently [7] for 2-connected triangulations, but no ones knows at the moment 
whether it will remain isolated or is just the tip of another iceberg. 

3.6. Algebraic series: some questions. We begin with three simple 
classes of objects that have an algebraic GF, but for reasons that remain mysteri- 
ous. We then discuss a possible criterion (or necessary condition) for N- algebraicity, 
and finally the algebraicity of certain hypergeometric series. 
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3.6.1. Kreweras' words and walks on the quarter plane. Let Z be the set 

of words u on the alphabet {a, 6, c} such that for every prefix v of m, \v\a > \v\h 
and \v\a > \v\c- These words encode certain walks on the plane: these walks start 
at (0, 0), are made of three types of steps, a = (1, 1), b = (—1, 0) and c = (0,-1), 
and never leave the first quadrant of the plane, defined by a;, y > 0. The pumping 
lemma \52\ Thm. 4.7], applied to the word a"6"c", shows that the language C is 
not context-free. However, its generating function is algebraic. More precisely, let 
us denote by iij (n) the number of words m of £ of length n such that \u\a — \u\b = i 
and \u\a — \u\c = j. They correspond to walks of length n ending at position 
Then the associated three- variable generating function is 

^ — ' — t(u + u + u^w^l uvt 

where u = 1/u, v = 1/v and W = W{t) is the unique power series in t satisfying 
W = t{2 + W^). Moreover, the number of walks ending at (i,0) is remarkably 
simple: 

^ , 4"(2^ + l) /2A/3" + 2^^ 



(n + i + l)(2?i + 2i + 1) V i 

The latter formula was proved in 1965 by Kreweras, in a fairly comphcated way [57] , 
This rather mysterious result has attracted the attention of several combinatori- 
ahsts since its pubhcation [HI [48l [64] . The first combinatorial explanation of the 
above formula (in the case j = 0) has just been found by Bernardi [7]. 

Walks in the quarter plane do not always have an algebraic GF: for instance, 
the number of square lattice walks (with North, South, East and West steps) of 
size 2n that start and end at (0, 0) and remain in the quarter plane is 

1 f2n + 2\ ^ 42"+ 1 



(2n-|- l)(2n-f 4) \ n + l 

and this asymptotic behaviour prevents the corresponding generating function from 
being algebraic (see l|16p ). The above formula is easily proved by looking at the 
projections of the walk onto the horizontal and vertical axes. 

3.6.2. Walks on the slit plane. Take now any finite set of steps 5 C Z x 
{ — 1,0,1} (we say that these steps have small height variations). Let Sij{n) be 
the number of walks of length n that start from the origin, consist of steps of S, 
never return to the non-positive horizontal axis {(— fc, 0), k > 0}, and end at (i, j). 
Let S{u, v; t) be the associated generating function: 

S{u,v]t) — Sij{n)u^v-'t"' . 

z j(zZ,n>0 

Then this series is always algebraic, as well as the series Sij{t) X]ri^ij('^)^" 
that counts walks ending at OEQ]. For instance, when S is formed of the 
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usual square lattice steps (North, South, West and East), then 

(1 - 2t(l + u) + V]~4t)^^^ (1 + 2t{l -u) + y/miy^'^ 



S{u^ v; t) 



1 — t{u + u + V + v) 



with u ~ \/u and v = l/v. Moreover, the number of walks ending at certain 
specific points is remarkably simple. For instance: 



si,o(2n + 1) = C; 



2ri+l, 



so,i(2n + 1) = 4"C„, s-i,i(2n) = C2. 



where C„ = (^^) / (n+1) is the nth Catalan number, which counts binary trees (fT7| . 
Dyck words, and numerous other combinatorial objects [ZH Ch. 6]. The first of 
these three identities has been proved combinatorially [4]. The others still defeat 
our understanding. 

3.6.3. Embedded binary trees. We consider again the complete binary trees 
met at the beginning of Section 13.21 Let us associate with each (inner) node of 
such a tree a label, equal to the difference between the number of right steps and 
the number of left steps one takes when going from the root to the node. In other 
words, the label of the node is its abscissa in the natural integer embedding of the 
tree (Fig.IEl). 





Figure 12. The integer embedding of a binary tree. 

Let Sj = Sj (i, u) be the generating function of binary trees counted by the number 
of nodes (variable t) and the number of nodes at abscissa j (variable u) . Then for all 
j G Z, this series is algebraic of degree (at most) 8 (while Sj{t, 1) is quadratic) |12j . 
Moreover, for j > 0, 



where 



r = i + tT\ 



z = t- 



il + Z' 



{i-z + z-^Y 

and yU = ^{t,u) is the unique formal power series in t satisfying 

Z{l + tiZ)\l + nZ^){l + liZ'') 



= (u- 1) 



(1 + z)2(i + z + z^){i- zY{i - n^z^) ■ 
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Why is that so? This algebraicity property holds as well for other families of 
labelled trees |121I24). From these series, one can derive certain limit results on the 
distribution of the number of nodes at abscissa [An^/^J in a random tree with n 
nodes [12] . These results provide some information about the law of the integrated 
super- Brownian excursion \12\ I21|. 



3.6.4. N-algebraicity. N-algebraic series have been defined in Section 13.21 in 
terms of positive proper algebraic systems. The author has been unable to find in 
the literature a criterion, or even a necessary condition for an algebraic series with 
coefficients in N to be N-algebraic. Nor even an algebraic series with coefficients 
in N that would not be N-algebraic (together with a proof of this statement...). 

A partial answer could be provided by the study of the possible asymptotic be- 
haviour of coefficients of N-algebraic series. It is very likely that not all behaviours 
of the form (fT6|) are possible. An important result in this direction states that, if 
a proper positive system lfT8|) is strongly connected, the nth coefficient of, say, 
follows the general pattern (US]), but with d = -3/2 |43l Thm. VII.7]. The system 
is strongly connected if, roughly speaking, the expression of every series Ai involves 
(possibly after a few iterations of the system) every other series Aj . For instance, 
the system defining the walks ending at in Section [3.4. II reads 

Mo = t^{l + Mof and Wo^Mo{2 + Wo). 

This system is not strongly connected, as AIq does not involve Wq. Accordingly, 
the number of 2n-step walks returning to is (^") ~ K4"n~^/^. 

If one can rule out the possibility that d ~ —5/2 for N-algebraic series, then 
this will prove that most map generating functions are not N-algebraic (see the 
examples in Section! 



3.6.5. Some algebraic hypergeometric series. Consider the following series: 

where oi, . . . , a^, 6i, . . . , 6e are positive integers. This series is algebraic for some 
values of the a^'s and 6j's, as shown by the case 



n>0 



Can we describe all algebraic cases? Well, one can easily obtain some necessary 
conditions on the sequences a and b by looking at the asymptotics of /„. First, 
an algebraic power series has a finite, positive radius of convergence (unless it is a 
polynomial). This, combined with Stirling's formula, gives at once 

ai H \-ad = bi-\ h 6e. (23) 

Moreover, by looking at the dominant term in the asymptotic behaviour of /„, and 
comparing with fT6|) , one obtains that either e = d, ore = d+ l. The case d = e 
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only gives the trivial solution F(t) = 1/(1 — i), and the complete answer to this 
problem is as follows |1H [69] : 

Theorem 3.11. Assume ((23]) holds, and F{t) ^ 1/(1 - t). The series F{t) is 
algebraic if and only if fn G N for all n and e = d + 1. 

Here are some algebraic instances: 

(6n)!(?i)! (10n)!(n)! (20n)!(n)! 

~ (3?i)!(2n)!2' ^ (5n)!(4n)!(2n)!' ~ (10n)!(7n)!(4?i)! ' 

The degree of these series is rather big: 12 [resp. 30] for the first [second] series 
above. This theorem provides a collection of algebraic series with nice integer 
coefficients: are these series N-algebraic? Do they count some interesting objects? 

Acknowledgements. The parts of this survey that do not deal exactly with the 
enumeration of combinatorial objects have often be influenced by discussions with 
some of my colleagues, including Frederique Bassino, Henri Cohen, PhiHppe Flajo- 
let, Frangois Loeser, Geraud Senizergues. Still, they should not be hold responsible 
for any of the flaws of this paper... 
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